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ABSTRACT 

A new chemical model is presented for the carbon-rich circumstellar envelope of 
the AGB star IRC-l-10216. The model includes shells of matter with densities that 
are enhanced relative to the surrounding circumstellar medium. The chemical model 
uses an updated reaction network including reactions from the RATE06 database and 
a more detailed anion chemistry. In particular, new mechanisms are considered for 
the formation of CN~, CsN" and C2H~, and for the reactions of hydrocarbon anions 
with atomic nitrogen and with the most abundant cations in the circumstellar envelope. 
New reactions involving are included which result in the production of significant 
amounts of C2H~ and CN~ in the inner envelope. The calculated radial molecular 
abundance profiles for the hydrocarbons C2H, C4H and CgH and the cyanopolyynes 
HC3N and HC5N show narrow peaks which are in better agreement with observations 
than previous models. Thus, the narrow rings observed in molecular microwave emission 
surrounding IRC-l-10216 are interpreted as arising in regions of the envelope where the 
gas and dust densities are greater than the surrounding circumstellar medium. Our 
models show that CN~ and C2H~ may be detectable in IRC+10216 despite the very 
low theorised radiative electron attachment rates of their parent neutral species. We 
also show that magnesium isocyanide (MgNC) can be formed in the outer envelope 
through radiative association involving Mg+ and the cyanopolyyne species. 

Subject headings: astrochemistry - circumstellar matter - stars: individual: IRC+10216 - ISM: 
molecules 

1. Introduction detected in the expanding circumstellar enve- 

lope (CSE), which extends more than 10,000 
The carbon-rich AGB star IRC+10216 is from the star. The outflow from the cen- 

among the richest known sources of molecules ^^^^ ^^^^ carbon-rich (with [C] / [O] > 1 ) , and 

in the sky. The star is nearing the end of its contains abundant stable molecules such as 

life and is rapidly losing mass to the inter- C2H2, HCN and CH4, which are formed in 

stehar medium. Over 60 species have been ^j^^ j^^^ -^^^^ ^^^-^^^ ^^^^ g^^H^^ g^^f^^^ 
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dLafont et alJll982l ). As the density falls with 
radius, interstellar UV photons penetrate the 
CSE and cause the photolysis of these 'par- 
ent' molecules. This results in the production 
of reactive species that subsequently partici- 
pate in a complex gas-phase chemistry to pro- 
duce new 'daughter' molecules in the outer 
envelope. Chemical models of this environ- 
ment (see iMillar et all bood . for example), 
are able to reproduce with good accuracy the 
observed column densities of outer-envelope 
species such as the carbon chains CnR (for 
n = 2,4,6,8). However, it has been sug- 
gested that such chemical models are incon- 
sistent with observed Cn H rotational e i nissiq n 
maps (see for example Guelin et al. 19991 ). 



The observations show C„H emission max- 
ima at around the same radius (for different 
n), whereas the results of chemical models 
show spatially separated peak abundances. 
The modeled molecular radial distributions 
are also broa der than observed. Using op- 
tical imaging, iMauron &: Huggin i (|200d l de- 
tected multiple dust shells in the envelope 
of IRC+10216, spaced at regular intervals 
of ~ 5" to 20". These shells were inter- 
preted as arising as a result of modulation 
of the stellar mass-loss rate, perhaps due 
to t he presence of a binar y stellar compan- 



Brown k Millail ^200i ) included density- 



ion. 

enhanced circumstellar dust shells in their 
chemical model for IRC-l-10216 in an attempt 
to make the model more physically realis- 
tic and to address the discrepancies between 
observed molecular emission maps and the 
results of previous chemical models. The ad- 
dition of dust shells modifies the radiation 
field inside the CSE and results in modeled 
molecular abundance profiles that are in bet- 
ter agreement with the observations. 



Fig ure [J shows the effect of the lBrown &: Millar 
(|2003l i density-enhanced dust shells on the 
modeled C2H, C4H and CqH radial pro- 
files in IRC-l-10216. In the presence of dust 



shells, the molecular abundance peaks move 
outwards and slightly closer together. The 
C2H profile FWHM narrows from 7" to 6" 
and moves f r om 8 " to 18". According to 
Guelin et al.l (|l999l l. the molecular emission 



maps from C2H, C4II and CgH all peak in a 
narrow circumstellar ring about 2" wide at a 
radius about 15" from the star. It is therefore 
clear that alt hough the radial distrib utions of 
the species in iBrown fc Millai] (|2003l Vs model 
represent a better match than previous mod- 
els, they are still too broad to fully explain 
the observed maps. In addition, CgH peaks 
at a larger radius than C2H and C4H, which 
is inconsistent with the observations. 

Dinh-V-Trung Liml ( 20081 ) mapped the 
molecular shells of HC3N and IIC5N in 
IRC-l-10216 at unprecedented angular reso- 
lution. The shells were found to be clumpy, 
co-spatial with each other and with a struc- 
ture closely matching the distribution of 
dust shells observed by iMauron &: Huggins 
(|200d ). ..t was concluded that the circum- 
stellar cyanopolyyne gas density distribu- 
tion matches that of the dust and therefore 
that the gas and dust are coupled. Shell- 
like stru ctures ^ 2" wide are presen t in the 
maps of Dinh-V-Trung Sz Lim ( 20081 ). which 
are inconsistent with the ~ 10" wide abun- 
dance profiles c alcula ted in the model of 
Brown fc Millail hooA ). 

In this article, a new chemical model for 
IRC-l-10216 is presented which builds on pre- 
vious models and includes density-enhanced 
sh ells of gas in ad d ition to the dust shells 
of iBrown fc Millail hooA ). Thus, the new 
model assumes kinematical coupling between 
the gas and dust in the CSE. The aim is 
to model the chemistry of the narrow shell- 
like struct ures observed in the I RC-l- 10216 
outflow by binh-V-Trung Liml too4 ) and 



Guelin et al.l (|l999l ) to test the idea that cou- 
pling between the gas and dust plays an im- 
portant role in determining the morphology of 



2 



E 0.15 




20 

Radius (arcsec) 









A C.H 




A 







Radius (arcsec) 



Fig. 1. — Radial molecular abundance profiles 
for C2H, C4H and CeH in IRC+10216, cal- 
culated in a model with no density-enhanced 
shells (top), a nd with the den si ty-en hanced 
dust shells of brown fc ^200^ (bot- 

tom). The distance to IRC+10216 has been 
assumed to be 130 pc. Details of the chemi- 
cal model used to calculate these profiles are 
given in Section [2j 



the molecular distributions. 

The molecular anions C4H , CgH , CgH 
and CsN" have been recently detected in 
the envelope of IRC+10216 (for a review of 
the astroph ysical anion detect i ons r eported 



far, see Herbst fc Osamura 20081 ). Fol 



so 

lowing these discoveries we also expand upon 
the anion chem i stry i n IRC +10216 studied 
by iMillar et al.1 (j200nl . 120071 ). In particular, 
the C„N~ chemistry is reconsidered using the 
reaction rates between N-atoms and carbon- 



chain anions published by lEichelberger et al 
(l2007l ). The possibility that magnesium iso- 



cyanide (MgNC) is produced in the outer CSE 
by gas-phase chemistry is also examined. 

2. The chemical model 
2.1. Physics 

The new chemical m odel for IRC+1021 6 
is based on the model of lMillar etaP (|200d '). 
The underlying density distribution in the 
CSE is derived using a mass-loss rate of 
1.6x10"^ Mgyr-i (after iMen'shchikov et~ar 



2OOII ) , and assuming that the matter expands 
in a spherically symmetric outflow with a ve- 
locity of 14 km s~^. The resulting gas num- 
ber density profile n(r) falls as 1/r^, added to 
which are a series of step-like density enhance- 
ments of the form (3n{r). Based on the dust 
shell parameters deduced from scatter e d ligh t 
observations by Mauron &: Huggins (2000), 
each shell is 2" thick with an inter-shell spac- 
ing of 12", which corresponds to a timescale 
of approximately 530 years between successive 
episodes of enhanced mass-loss. The param- 
eter (3 is set to 5 for all shells in the present 
model (see Section [23]) • To convert between 
physical length units and angular distances on 
the sky, the dista nce to IRC+10216 is take n 
to be 130 pc (after IMen'shchikov etaPboOll ). 
The chosen shell parameters are not intended 
to provide an accurate representation of the 
dust shell structure observed in IRC+10216 
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(for which a 3-dimensional, time-dependent 
model of the gas and dust in the envelope 
would be required), but to permit the study 
of the general effects on the chemical model of 
the addition of density-enhanced gas and dust 
shells. The radial H2 gas density distribution 
used in the model is shown in Figure [5J The 
H2 is assumed to be completely self-shielded 
in the regions of interest in the CSE so that 
nn^i^) ~ n(r). 

The adopted temperature profile is based 
on an empirical fit to tVip gas Vinptir tpm- 
neiature profile derived by ICrosas &: Menten 
(119971) ■ and takes the form 



r(r) 



2.81 X 10^5 



3.34 X 10^ 



17 \ 1-05 



, 10 



(1) 

with a lower limit of 10 K to prevent the tem- 
perature becoming unrealistically low in the 
outer envelope. Parent species (with abun- 
dances shown in Tabled]) are injected into the 
model at the inner radius of = 10^^ cm, 
where the gas density is 8.6 x 10'' cm~^. At 
this radius the gas kinetic temperature is 575 
K and the interstellar radiation field is atten- 
uated by an effective radial extinction of 47 
magnitudes in the V band. Initia l abun dances 



are the same as in Millar et al 



(|2000l l. with 
the addition of Mg (see Section [3T3]) . and with 
C2H2 and HCN abundances taken from the 



recent measurements by Fonfria et al. ( 20081 ). 
The standard interstellar radiation field 
is assumed to impinge on the 



(|Drainelll97; 

outside of the circumstellar envelope from all 
directions. The extinction is calculated for 
the underlying 1/r^ density distribution of 
the C SE using the approach of lJura &: Morris 
(|l98lh . Additional terms are added to the 
extinction due to the contributions of the 
density-enhanced shells (for further details 
of the extinction calculation in the pre sence 
of dust shells, see lSrown fc MillaillioO^ ). 



2.2. Chemistry 

The chemical reaction 
on that used by 



Petrie et al 



networ k is based 
(I2OO3I I and 



Millar et al.l (120071 ). The reaction rates have 
been updated to be consistent with those 
i n the (dipole-enhan ced) RATE06 database 



(|Woodall et al 
actions from RATE06 



20071 ). 



Additional new re- 
for those species in 
the IPetrie et al.l (|2003l ) model have also been 
added. The following reactions were deleted 
from the reaction network in order to increase 
computational speed, with negligible effect on 
the chemistry of the species of interest in this 
study: those with activation energies greater 
than 300 K, those reactions involving HaO^ 
as a reagent (except for HsO"*" + e~), those 
with as a reagent (except for + e~ 
and + H2), and any reactions involving 
C", S-, C0+, HOC+, C2H+ and C2H5. 

Carbon-chain species, hydrocarbons, cyanopolyynes 
and their associated anions and cations are 
of principal importance for the chemistry of 
this study. The following species are among 
those included in the chemical model: carbon 
chains C^ (n = 1 - 23), C+ (n = 1 - 23), 



n 



n 



HC2n-lS , 

= 2 - 23), 
= 1-4), 
= 1-3), 



3 - 23); C^t^iS (n = 1 - 3); 



1 — 3); hydrocarbons C„H(+/ 



(+) 



1 



3 

CnH^ 



1 



23), C„H3 
23), C„H4 
9), C,H+ 

) 



(n 
(n 

(n = 1 

- 9); and cyanopolyynes C2n-iN("^/ 
1 - 11), HC2„_iN (n = 1 - 12), 
HC2n-iN+ (n = 1 - 11) H2C2n_iN+, (n = 
1-11), H3C2„-iN+ (n = 2-5). 

Due to the current interest in molecular 



anions in IRC+10216 (e.o.lMillar et al 



Remiian et al.l 120071: iThaddeus et al 



Cordiner et al. 



200. 



2007; 



2008 



the anion chemistry 
has been extended to include C^H" down to 



n 



2, utilizing t he radiative electron , attach- 
ment rates from Herbst fc Osamura ( 20081 ). 
The following additional CN~ and C2H~ for- 
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Table 1 

Initial fractional abundances of parent species relative to H2 



Species 


Initial abundance 


He 


1.5x10^1 


C2H2 


8.0x10"^ 


CH4 


2.0x10-6 


H2S 


1.0x10-6 


HON 


2.0x10-^ 


NH3 


2.0x10-6 


CO 


6.0x10-"^ 


cs 


4.0x10-6 


N2 


2.0x10-"^ 


Mg 


1.0x10-^ 



mation reactions have been included: 
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Fig. 2. — Radial H2 density profile used in the 
model. 



HCN + 



C2H2 + H 



CN- + H2 



C2II + II2 



(2) 
(3) 



The rate coefficient used for Equation ^ 
(3.8 x 10- ^ cm^ s~^), is an estimate taken from 



Prasad HuntressI (|l98d l. For Equat: 
the rate coefficient (4.42 x 10 



ion 



Q 'H — 1 \ 

cm s I was 



measu red experimentally by iMackav et al 
(|l977l ). These proton-transfer reactions are 



likely to be important in the model due to 
the large HCN and C2H2 abundances in the 
stellar outflow. H- is produced in the model 
mainly by the cosmic-ray (CR) dissociation of 
H2 (H2 + CR — > H+ + H-). This reaction 
is slow ( with a rate coefficient of 3 .9 x 10-^^ 



cm^ s 



1. 



Prasad fc HuntressI Eosol 'l. but pro- 



vides the main source of H- in the inner enve- 
lope. Interior to the second density-enhanced 
shell (r < 2 X 10^6 ^m), the modeled H- 
abundance is about 10^^ cm-^. Other re- 
actions similar__toEquati^ ([3]) were stud- 
ied by Mackav et al. ( 19771 ). who found that 
many different molecular anions could recieve 
a proton from C2H2 at rapid rates (~ 10-^ 
cm^ 3-^), resulting in the production of C2H-. 
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However, it is presently unknown whether re- 
actions occur between C2H2 and the carbon 
chain anions C~ and C^H". In hght of an 
observational upper Hmit f o r A^( C2H~) in 



IRC+10216, ICordiner et al.1 (|2008l 1 deduced 



that these reactions probably do not proceed 
rapidly. 

Dominant anion destruction mechanisms 
are by reaction with H and and by 
photodetachment. Photodetachment rates 
were calcula t ed ac cording to Equation (2) of 
Millar et aD (|2007l ). As a result of its large 
electron detachment energy, the CN~ pho- 
todetachment rate thus calculated is ~ 100 
times less than the value in the RATE06 
database. 

Carbon chain anions (n = 2 — 7), and 
CnH~ (n = 2, 4, 6) have been shown to react 
with atomic nitrogen and result in the forma- 
tion of products tha t include the nitrile anion s 
C„N- (n = 1,3,5) dEichelberger et al.ll2007h . 
We have included these reactions in the model 
as part of our investigation into the possible 
mechanisms for the formation of CN~ and the 
recen tly discovered CsN" (IThaddeus et al 



20081) . Branching r a tios w ere calculated from 
Eichelberger et al.l (|2007l Vs original experi- 
mental data by Veronica Bierbaum (private 
communication). The branching ratio for the 
associative electron detachment (AED) prod- 
uct channel could not be derived from the 
experimental data, so we have assumed, arbi- 
trarily, a ratio of 0.5. This assumption consti- 
tutes potentially the most significant source 
of error in these reaction rates. 

Mutual neutralisation reactions between 



anions and cations were shown by lLepp &: Dalgarno 
(|l988l ) ;o have important effects on interstel- 
lar chemistry. Thus, reactions of the kind 



X + Y 



(4) 



have been included for the twenty most 
abundant cations (X"^), and for all anions 
(Y~), in the model, with a rate coeffi- 



cient of 7.5x 10~^(r/ 300)~°-^ cm^s"^ (see 
Harada Sz Herbst 20081 . for example). 

The final chemical network contains 426 
gas-phase species coupled by 5539 reactions. 

2.3. Evolution of the shells 

Initially the rate equations are solved as 
a function of rad ius (in the same way as 
Millar et all boool ). starting from = 10^^ 
cm and moving out to = 3 x 10^^ cm where 
photodissociation destroys all the molecules 
(apart from H2). Then in a separate routine 
the chemical abundances inside the density- 
enhanced shells are calculated as a function 
of radius. The chemical rate equations of 
a density-enhanced packet lying halfway be- 
tween a shell's inner and outer radius (at 
a radius rp with density {(3 + l)n(rp)), are 
solved, starting from and moving out to 
2 X 10^^ cm where the density is sufficiently 
low that the shell no longer makes any signifi- 
cant contribution to the total amount of mat- 
ter in the model. The motion of the density- 
enhanced shells in the outflow are followed 
so that they move outward over time, syn- 
chronised with the outward radial motion of 
the dense packet. The radiation field is re- 
calculated for the dense packet at each time- 
step, taking into account the new positions 
of all the shells. The density factor (/?), the 
inter-shell spacing and the shell thickness are 
identical for every shell so that the chemical 
abundances calculated for the dense packet 
at a given radius represent the abundances 
in a density-enhanced shell centered at that 
radius. 

2.4. Model variations 

The effects on the model results of vari- 
ations in the mass-loss rate, radiation field 
strength, gas-to-dust ratio, shell thickness, 
inter-shell spacing, density-enhancement fac- 
tor and stellar distance have been analysed. 
The results presented in this study are for the 
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model that uses the parameters that we be- 
heve best match the observational constraints. 
Modification of these parameters, particularly 
those that affect the radiation field strength, 
can significantly alter the radial abundance 
profiles calculated by the model. However, 
under such circumstances the main conclu- 
sions of this study remain the same. The lo- 
cation of the 15" density-enhanced shell has 
been fixed in order to best reproduce the ob- 
servational data, which may be considered a 
contrivai ice of the model. Howeve r, in Fig- 
ure 8 of iMauron &: Huggin this can 
be identified as the radius at which the first 
distinct dust shells occur. 

To permit comparison between this model 
and previous models published in the litera- 
ture (which do not include density-enhanced 
gas and dust shells), we have also run the 
model without any density-enhanced shells 
{i.e. /3 = 0) and also wi th the density- 
enhanced dust shells of (Brown &: Millar 



20031 ) to produce the data shown in Figure 



m 



2.5. Molecular excitation 

To facilitate comparison of the model re- 
sults with observed maps of molecular mi- 
crowave emission line flux, the rotational ex- 
citation of some of the molecules of inter- 
est has been calculated function of ra- 
dius using a modified versi on of the mmline 



comp uter code (described bv lJusttanont et al 
19941 ). The central stellar radius was taken 
to be 3.1x10^^ cm and the temperatu r e 265 



K (derived from Men'shchikov et aD I2OO1I ). 
T he dust opacity was taken from Figure 6 



of iMen'shchikov et al.1 (|200ll ). Due to the 



lack of published collisional excitation rates 
and vibrational transition strengths, only a 
rough estimate of the rotational excitation is 
possible for most molecule s. For C2H and 



HC3N rates from Green Sz Chapman ( 19781 ) 
are used. For CgH and CgH" the rates have 
been extrapolated up to J = 31. The colli- 
sional rates used are for closed electronic-shell 
species so the spectroscopic structure of the 
corresponding (closed electronic-shell) anions 
has been used for the calculations of the exci- 
tation of the open-shell C2H, C4H and CeH 
radicals. This approximation is reasonable 
because the structure of these hydrocarbon 
anions and neutrals are very similar. Only 
the ^Hi/2 states have been considered for the 
neutral hydrocarbons; the population of the 
^n3/2 states are not expected to significantly 
affect the relative populations of the states of 
interest here. Rotational Einstein A coeffi- 
cient s were calcul ated using dipole-moment s 
from lWoon I (|l995l ^ and lBlanksbv et al.l (|200lh 
for the neutrals and the anions, respectively. 
IR pumping of rotational levels has been cal- 
culated through consideration of the radiative 
excitation of a single vibration al state ^ 1 



jum a bove the ground state (c/. iBieging et al 



1984 1. C2H has been calculated to have a 



strong vib rational transition (A = 0.6 s~^) at 
12.5 iJ,m (jTarroni Sz Carteii 12004 ). which we 
assume to also occur in C2H~. The vibra- 
tional spectra of C4H, CgH and their associ- 
ated anions are less well known. For these 
species a transition has been assumed to oc- 
cur at 12.5 ^m with A = 1 s~^. IR pumping 
has a significant effect on the molecular exci- 
tation, but changes in the vibrational transi- 
tion wavelengths and Einstein A coefficients 
by up to an order of magnitude do not signif- 
icantly affect the results of the present study. 
To assess the impact of errors in the collisional 
excitation rates on the calculated molecular 
emission profiles, the rates were varied by an 
order of magnitude either way. The overall 
features of the emission profiles remained the 
same. 



C2H- the HCO+ rates from lFlowerl (|l999l l are 

used, and for C4H, C4H~, CgH and CgH" the 
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3. Results 



3.1. Molecular radial abundance and 
intensity distributions 

The upper panel of Figure [3] shows the cal- 
culated radial abundance profiles for the hy- 
drocarbons C^H (n = 2,4,6). Comparison 
with Figure [1] shows the impact of the density- 
enhanced shells on these species. A prominent 
effect of the shells is to raise of the abundances 
relative to the surrounding CSE due to the in- 
creased gas density. The photon-induced hy- 
drocarbon chemistry is suppressed very near 
to the star by the dust shell at a radius of 
r = 1" which provides additional shielding 
of C2H2 from photodissociation. The abun- 
dances all reach a peak in the shell at r = 15". 
However, this is not where CeH reaches its 
greatest abundance, which occurs in the third 
density-enhanced shell at r = 29". As shown 
in the middle panel of Figure [U HC3N and 
HC5N both reach their greatest abundances in 
the 15" shell. Cosmic-ray-induced chemistry 
also results in the synthesis of a significant 
amount of HC3N in the innermost (r = 1") 
sheU (through the reaction HCN + CRPHOT 
— ^ CN H, followed by CN + C2H2 — > 
HC3N + H). The anions C4H-, CeH" and 
CsH" (shown in Figure E]), do not reach their 
maximum abundances in the 15" shell; C4H~ 
peaks at a similar radius to the model with 
no shells whereas CgH" and CgH" peak in 
the 29" sheU. 

Molecular abundances tend to peak within 
the density-enhanced shells because the in- 
creased density raises the abundances of 
chemical reagents which drives the chemistry 
at a faster rate. It is not always the case, 
however, that this raises the abundances of 
daughter species, as can be seen in Figure [3] 
where the C2H abundance is reduced in the 
r = 29" shell due to the increased densities 
of atoms and ions (including C, C^ and N), 
that it reacts with. 




10 20 30 40 

Radius (arcsec) 

Fig. 3. — C2H, C4H and CgH modeled radial 
abundance profiles (top) and 3 mm emission 
intensity profiles (multiplied by constant scal- 
ing factors for display) (bottom). 
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Fig. 4. — Modeled radial abundance pro- 
files for HC3N and HC5N in the absence 
of density-enhanced shells (top) and with 
density-enhanced shells (middle). The bot- 
tom panel shows the modeled HC3N J = 5 — 4 
intensity profile calculated in the presence of 
density-enhanced shells. 



A comparison of the abundance profiles 
in the models with and without density- 
enhanced shells shows that, in general, the 
increased shielding of the CSE from interstel- 
lar UV by the dust shells (which inhibits the 
photochemistry), causes photodissociation of 
parent species to be less efficient in the inner 
regions and causes the daughter abundances 
to rise more slowly with radius, moving the 
profile maxima outwards. 

The calculated emission intensities for 
C2H, C4H, CqH, CeH- and HC3N (shown 
in the lower panels of Figures [3l H] and [5l re- 
spectively), are greatest within the r = 15" 
density-enhanced shell. C4H~ reaches maxi- 
mum intensity near r = 8". The differences 
between the abundance profiles and the emis- 
sion profiles for CeH and CeH" are partic- 
ularly notable. Because the 3 mm emission 
from these species originates from a high ro- 
tational level (around J = 30) , the strength of 
the emission is highly dependent on the rate of 
collisional excitation, and therefore the den- 
sity. Thus, inside the density-enhanced shells 
the lower J levels tend to become depopulated 
in favour of the higher levels. 

For C2H, C4H, CeH and their correspond- 
ing anions, two-dimensional emission maps 
have been constructed from the radial emis- 
sion profiles by integration along (pencil- 
beam) lines of sight through the CSE. This 
calculation assumes spherical symmetry of the 
CSE and therefore cannot show any of the az- 
imuthal structure in observed emission maps. 
Nevertheless, these maps (shown in Figure 
[6]), provide a useful means for comparing the 
main features in the modeled and observed 
molecular emission. 

C2H and CeH show a strong, narrow emis- 
sion ring centered on 15", in ex cellent agree- 



Guelin et al 



ment with observed maps of 

(|l999l ). [n addition to the observed 15" ring, 
the modeled C4H emission map shows a thick, 
strong ring at around r = 8'. The patterns for 



9 




0.0009 
0.0008 
0.0007 
0.0006 
0.0005 
0.0004 
0.0003 
0.0002 
0.0001 




Radius (arcsec) 



C.H" 




Radius (arcsec) 



Fig. 5. — Modeled radial abundance pro- 
files for the hydrocarbon anions C4H , CgH 
and CgH" in the absence of density-enhanced 
shells (top) and with density-enhanced shells 
(middle). C4H~ and CeH" 3 mm emission in- 
tensity profiles are shown in the bottom panel. 



C4H~ and CgH" are similar to their parent 
neutrals whereas C2H~ is markedly different: 
it has a centrally-peaked emission map be- 
cause it is produced predominantly by the re- 
action of H~ with C2H2 — the abundances of 
which are greatest in the (well shielded) inner 
CSE — rather than by electron attachment 
to its parent neutral. 

3.2. New anion results 

Column densities for all anions included in 
the chemical model are given in Table [2j 

Cor diner et al. ( 20081 ) identified that the 
reaction of H" with HCN (Equation ([2])) 
could dominate the synthesis of CN~ in 
IRC-l-10216. However, the present study 
shows that reactions between atomic nitro- 
gen and the carbon chain anions may provide 
a greater source of CN~. Upon inclusion of 
the reactions N -|- C~ and N -|- CmH~ into 
the chemical network (see Section 12. 2p . the 
CN~ column density is raised from 1.3 x 10^'^ 
to 1.3 X 10^2 cm~2. Due to the fact that Cf 
is the most abundant anion involved in these 
reactions, the dominant formation reaction 
for CN- is N + Cf — > CN- + Ce. These 
results are dependent on the values of the 
branching ratios assumed for the associative 
electron detachment reactions. Even if the 
AED branching ratio is as large as 0.9 (in- 
stead of the assumed value of 0.5), the CN~ 
column density is calculated to be 2.7 x 10^^ 
cm~^ and N -|- Cj is still the dominant reac- 
tion for CN~ production. 

CsN" was detected for the first time in 



IRC-M0216 by iThaddeus et~aD ^2Q04 ). who 
report a column density of 1.6 x 10^^ cm~^. 
If CsN" is assumed to be formed only by ra- 
diative electron attachment to C3N, the mod- 
eled column density is 1.1 x 10^'^ cm~^, with a 
corresponding anion-to-neutral ratio of 0.02%. 
This would imply that the electron attach- 
ment rate used (calculated by lPetrie &: Herbst 
19971 ). is too small. However, upon inclu- 
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Table 2 

Calculated column densities of negatively- charged species 



Species N (cm ^) Species A'' (cm ^) Species A'' (cm ^) 



C3 


3.9cl0 


C2H- 


5.5el0 


CN- 


1.3cl2 


C4" 


7.5el2 


C3H- 


1.4e9 


C3N- 


g.Oell 


C5- 


3.2el3 


C4H- 


2.5el3 


C5N- 


1.7el3 


Ce 


5.3el3 


C5H- 


1.2el3 


C7N- 


7.2el2 


C7 


1.6el4 


CeH- 


9.6el3 


CgN- 


3.6ell 


Cs 


9.6el3 


C7H- 


5.2el3 


CiiN- 


3.7el0 




1 1^1/1 

i.iei4 




z.Ueio 


Ul3l\ 


z.oeiu 


•^10 


4.0el3 


CgH- 


1.6el3 


C15N- 


1.9el0 


cri 


3.1el3 


CioH- 


5.2el2 


C17N- 


1.3el0 


^12 


8.3el2 


CiiH- 


1.4el3 


C19N- 


8.5e9 


^13 


1.5el3 


C12H- 


3.1el2 


C21N- 


5.3e9 


Ci4 


3.9el2 


C13H- 


8.4el2 






crs 


4.6el2 


C14H- 


1.4el2 


CH2CN- 


2.5e6 


^16 


1.8el2 


C15H- 


6.4el2 








2.5el2 


CieH" 


l.lel2 


H- 


3.3e8 


^18 


l.lel2 


C17H- 


4.6el2 






^19 


1.5el2 


CisH- 


9.4ell 


e~ 


4.8el5 


•^20 


6.9ell 


C19H- 


3.2el2 






*^21 


1.0el2 


C20H- 


7.8ell 






^22 


4.2ell 


C21H- 


2.1el2 






*^23 


S.lell 


C22H- 


5. Sell 








C23H- 


9.9el2 







sion of reactions between nitrogen atoms and 
anions, the modeled CsN" column density 
rises to 9.0 x 10^^ cm~^ (which corresponds 
to an anion-to-neutral ratio of 1.3%), in rea- 
sonable agreement with the observed value. 
The dominant reactions involved in the pro- 
duction of CsN" in our present model are N -|- 
C~ (n = 5,6,7). If these reactions do indeed 
dominate the CsN" synthesis then the rate of 
radiative electron attachment to C3N calcu- 



lated by iPetrie fc Herbsti (|l997l ) may still be 
accurate. 

Cyanopolyynes C„N (for n > 4), by 
analogy with the structurally similar lin- 
ear hydrocarbons CwH (n > 5) studied by 



( Herbst Qsamura 20081 ). have been 



as- 
sumed to undergo rapid radiative stabilisation 
upon attachment of a free electron so that the 
rate of radiative electron attachment used in 
the model is 1.25 x 10"'^(T/300)"°-^ cm^^ s"^ 
Consequently, CsN^ and CrN^ have large 



modeled column densities of 1.7 x 10^^ and 
7.2 X lO^'^ cm~^ respectively, which correspond 
to about 8% of the abundance of their neutral 
parents C5N and C7N. Radiative electron at- 
tachment dominates the production of these 
anions in the model, but it seems plausible 
that reactions between atomic nitrogen and 
carbon -chain anions longe r than those stud- 
ied by lEichelberger et al.l ()2007l ) would also 
result in fragmentation of the carbon chain. 
Such fragme ntation was observed in t he ex- 
periments by lEichelberger et al.l (|2007l ) (who 
studied carbon chain lengths up to only n = 7 
carbon atoms), which gave rise to products 
with carbon chains with lengths from 1 to 
n. If the fragmentation of carbon-chain an- 
ions with n > 7 was included in our chemical 
model, there would be a significant increase 
in the production rates of the nitrile anions 
C2n-iN^ (n = 1 — 4), so that these reactions 
could potentially dominate the production of 
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Fig. 6. — Emission maps for hydrocarbons 
and their anions, calculated from the mod- 
eled radial intensity profiles assuming spheri- 
cal symmetry of the CSE. The grey scale is set 
to one at the respective emission maxima of 
each plot. Spatial units are arcseconds from 
the central star. 



C5N" and C7N" as well. 

Dense shells in the inner envelope have a 
significant impact on the C2H~ radial abun- 
dance profile, as shown in Figure [71 C2H^ 
is concentrated in the region of the envelope 
around the r = 1" shell. The presence of this 
shell raises the C2H~ column density from 
2.3 X 10^° to 5.5 X 10^° cm-2. This anion 
is unusual in the model because its rate of 
formation by radiative electron attachment is 
very slow; it is produced predominantly in the 
inner envelope as a result of Equation ^ (H~ 
+ C2H2 — > C2H~ + H2). This may have im- 
portant observational consequences because 
the kinetic temperature of the gas is higher at 
such radii, which would cause the molecule to 
exist in a higher state of rotational excitation. 
Direct comparison of (unresolved) single-dish 
microwave observations of C2H~ and C2H 
may be made more difficult by the difference 
in the distributions and therefore the tele- 
scope beam-filling factors of these two species. 
Detailed microwave observations of the spa- 
tial distributions of C2H~ are clearly required 
in order to confirm this result and determine 
whether Equation ([3]) is indeed the dominant 
production mechanism for this species. 

In the new model for the circumstellar en- 
velope of IRC-l-10216, the respective anion-to- 
neutral column density ratios for C4H, CeH 
and CgH are calculated to be 1.4%, 7.4% and 
4.5%. The corres ponding observational r atios 
are 0. 024%, 6.3% (ICernicharo et al.l boO';]'! and 



26% (|R,emiian et al.l 120071 ). As highlighted 



above for C2H, obtaining accurate observa- 
tional ratios is hindered by the possibility that 
the anions and neutrals have different spa- 
tial distributions within the CSE. Neverthe- 
less, it is clear that our model still severely 
over-estimates the amount of C4H~ compare d 
with C4H (see also lHerbst fc: Osamurall2008l ). 
It may be the case that the C4H radiative elec- 
tron attachment rate calculated by Herbst is 
too large. Alternatively, C4H~ may be de- 
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stroyed more rapidly than our model currently 
predicts, for example, by reaction with HCN 
or C2H2 (of. Equations ([2]) and ([3])). In 
that case, such reactions would have to be 
much more rapid for C4H~ than for CgH^ 
and CgH", which could conceivably result 
from the different electron binding energies 
of these species (Herbst, private communi- 
cation). Further laboratory and/or theoret- 
ical studies are required in order to determine 
whether the organic anions considered here re- 
act with HCN, C2H2 or other possible proton- 
donating molecules. 

Anion-to-neutral ratios of the C^i and CnH. 
species are typically slightly greater in the 
present model compared to the no-shells (NS) 
model. This is largely attributable to the in- 
creased electron abundances inside the domi- 
nant 15" shell which leads to increased rates 
of electron attachment. Figure [8] shows the 
electron abundances function of radius. 
The density-enhanced shells cause increased 
shielding of the gas from photoionisation 
which results in electron abundances which 
are a factor ^ 2 lower between r = 20" and 
100". Inside the shells at these radii, contrary 
to the 15" shell, the electron densities are 
up to an order of magnitude lower than the 
surrounding CSE. Beyond ~ 100" where the 
density becomes very low the shells begin to 
have a negligible effect on the photoionisation 
rate as shown for atomic carbon in Figure [8j 
Accordingly, the electron abundances of the 
models with and without shells converge at 
this radius. 




10 20 

Radius (arcsec) 



Fig. 7. — Modeled C2ll~ number density 
profiles with and without density-enhanced 
shells. The C2H abundance is also shown, 
multiplied by 10~^ for display. 




3.3. MgNC 

A number of metal-containing molecules 
have been detected in IRC-M0216. While 
some of these, for example NaCl, KCl and 
AlCl, are expected to be abundant in the 
LTE region and have spatial distributions that 
peak on the star, other species, particularly 
the cyanides and isocyanides such as MgCN, 



Fig. 8. — Electron number density profiles 
for models with and without density-enhanced 
shells. The atomic carbon photoionisation 
rates (C -1- hi') for the two models are also 
plotted. 
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MgNC and AINC, have shell- like distributions 
on scales very sim ilar to the cyanopolyynes 
(jZiurvs et al.ll2002l ). Two possibilities exist to 
account for the presence of these species in 
the outer envelope: either they form in the 
gas phase from metals that have not been in- 
corporated into dust in the inner envelope or 
they have been released from grains by erosion 
processes in the outer envelope. 

In order to test the hypothesis that metal 
isocyanides are formed by gas-phase chem- 
istry we have investigated the synthesis of 

MgNC via radiative association re actions be- 

tween Mg"^ and the cyanopolyynes (jPunbar Petrie 
2OO2I ): 



outward and to concentrate the abundance to 
the peaks of the gas density, in this case with 
almost equal peaks in the 15" and 29" shells. 
The total MgNC column density is calculated 
to be 5.7x10^'^ cm^^, compare d to observed 
values of fO.93-5'1 x 10^^ cm"^ (ICuelin et al " 



199,4 Iffighberger fc Ziu^booi ). The domi- 
nant formation reaction is with HC7N which, 
despite its smaller abundance has the largest 
radiative association rate coefficient (Equa- 
tion (El)), 6.59xl0-9(r/300)- 



i). 



-0.47 



cm 



Mg+ + HCnN — > HCnNMg+ + hz^ (5) 

(with n = 3,5,7,9), followed by dissociative 
recombination reactions of the form: 



HCnNMg+ + e" 



MgNC + Cn_iH. (6) 



Figure shows the radial abundance (cm~^) 
of MgNC in both the presence and absence of 
shells. In these calculations we adopt an ini- 
tial Mg/H2 abundance ratio of 10~^ — the 
derived abundances and column densities are 
directly proportional to this ratio. The ef- 
fect of the density-enhanced shells is, once 
again, to move the peak of the distribution 



0.00018 ^ 
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0.00012 - 

0.0001 - 

8e-05 - 

66-05 - 

4e-05 - 

26-05 - 




■ MgNC with shells 

■ MgNC no shells 



L 




Fig. 9. — Modeled MgNC number density 
profiles with and without density-enhanced 
shells. 



(iDunbar fc Petrid 12002 

It was sugg ested bv bunbar Petriel (|2002l ) 
that the inclusion of Equation ([5]) into chemi- 
cal models for IRC+10216 might significantly 
reduce the abundances of the species HC„N 
and help to reconcile the discrepancy between 
modelled and observed HC„N column den- 
sities. In the present model, we find that 
the reaction of Mg+ with HC5N, HC7N and 
HCgN results in only ~ 1% reduction in their 
calculated column densities, which does not 
significantly improve the match with observa- 
tions. However, the depletion of these species 
is dependent (in a roughly linear fashion), on 
the initial Mg abundance employed. 

4. Discussion 

Compared to prev ious chemical n iodels 
of IRC-H0216 ( e.g. iNeiad fc MillaJ Il987|; 

Millar et al. 200C; 



Cherchneff et al. 



1993 



Brown fc Millail boosi : 



Agundez et al.l Imm . 



the model presented here is unique in the in- 
clusion of density-enhanced shells of gas and 
dust. The shells included in our model have 
physical par ameters similar to th e dust shells 
observed by Mauron Sz Huggins ( 200d ). As 
expected, the modeled column densities dif- 
fer from th o se ca. lculated in the model by 
Millar et al. ( 200d ) (referred to hereafter as 
MHB). Table [3] gives the column densities for 
these two models for 33 species for which ob- 
servational column densities have been pub- 
lished. To highlight the effect of the density- 
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Table 3 



Comparison between column densities from the present mo del, the present mod el 

INCLUDING NO DENSITY-ENHANCED SHELLS (NS), THE MODEL BY IMILLAR ET AL.I f 2000l ) 

(MHB), AND OBSERVATIONS 



Species 


Present 


NS 


MHB 


Obs. 


Rof. 


c 


2.5el6 


2.7el6 


1.0el6 


I.lel6 


1 




2.6el5 


2.6el5 


9.9el5 


7.9el4 


2 


C2H 


5.6el5 


6.3el5 


5.7el5 


3-5el5 


1 


CN 


2.2el5 


2.6el5 


l.OelS 


l.lel5 


2 


HCO+ 


3.1el2 


3 0el2 


2.4el2 


3el2 


I 


Cq 


3.4el4 


2.6el4 


6.5el4 


lel5 


I 


CqH 


1.6el4 


1.2el4 


1.4el4 


3-7el3 


1 3 


CqHo 


1.0el4 


5.6el3 


2,lel3 


2el3 


1 


CH2CN 


2.0el3 


1.7el3 


6.9el2 


8.4el2 


4 


C3H4 


1.9el2 


5.7ell 


4.0el2 


1.6el3 


4 


CH3CN 


3.9el2 


3.6el2 


3.4el2 


6-30el2 


1, 4 


C4H 


1.7el5 


1.8el5 


l.OelS 


2-9el5 


1 


C4H- 


2.5el3 


2.1el3 




7.1ell 


5 


MgNC 


5.7el3 


6.6el3 




0.93-5el3 


6, 7 


C4H2 


8.0el5 


6.1el5 


2.9el5 


3-20el2 


1 


C3N 


7.0el3 


5.7el3 


3.2el4 


2-4el4 


1 


C3N- 


9.0ell 


7.7ell 


4.0ell 


1.6el2 


8 


HC3N 


6.7el4 


4.8el4 


1.8el5 


l-2el5 


1 


CH2CHCN 


2.6ell 


6.6el0 


l.lell 


5.1el2 


4 


C2S 


4.1el3 


4.4el3 


3.5el3 


9-15ol3 


9, 10 


Cs 


5.1el4 


6.3el4 


7.5el4 


lol4 


1 


C5H 


1.9el4 


2.0el4 


8.7el3 


2-50el3 


1 


C3S 


1.5el4 


1.7el4 


6.7el2 


6-llel3 


9, 10 


CeH 


1.3el5 


1.5el5 


5.8el4 


6.6el3 


5 


CgH- 


9.6el3 


l.lel4 




4.1el2 


5 


C5N 


2.4el4 


2.5el4 


1.4el4 


3-6el2 


1, 3 


HC5N 


3.5el5 


2.7el5 


7.1el4 


2-3el4 


1 


C7H 


3.7el4 


3.2el4 


4.5el3 


l-2el2 


1, 3 


C5S 


1.8el4 


1.7el4 




2.5el3 


10 


CgH 


4.1el4 


4.9el4 


l.lel4 


5el2 


1 


CgH- 


2.0el3 


2.1el3 


2.7el3 


2el2 


11 


HC7N 


1.4el5 


1.3el5 


2.2el4 


lel4 


1 


HC9N 


8.5el4 


8.0el4 


5.8el3 


3cl3 


1 



Ref erence s. — (1) see refere nces i n Table 5 of [ M illar e t al.l 
1I2OOOI'): (21 iBakker etHI l|l997h : (31 ICernicharo et al.l teOOOl '): 
(4) lAgundez et al.l l l200a) : fS) ICernicharo et al.l ||2007I)" (6) 
Guelin et al.l lll995l): (7) iHighberger fc Ziurvd ||2003^ : (8) 
Thaddcus et al I ll200i): f 9TTCernicharo et al .1 (Il987l1 : (10) 
Bel l ct a l. (JjgJ); Vll nRemiian et al.l 1120071 1. 
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enhanced shells on the chemistry and to per- 
mit comparison with models without density- 
enhanced shells, the column densities from 
the present model with no density-enhanced 
shells (referred to as NS), are also given. The 
column densities calculated by the three dif- 
ferent models are generally in good agreement 
with observations, especially given the com- 
plex morphology of the source. 

Notable differences between the MHB and 
NS models include a substantial reduction 
in the HC3N/HC5N column density ratio in 
the new model due to the enhanced pho- 
todissociation rate used for HC3N (from the 
RATE06 database). The C3N abundance 
is reduced as well because it is produced 
predominantly from HC3N (via an alterna- 
tive photodissociation channel). The carbon 
chains CeH, C7H and CgH, the cyanopolyynes 
HC5N, HC7N and HC9N and also CN are 
significantly more abundant in the present 
models than MHB. This is primarily due to 
the increased initial abun dances of the parent 
species C2H2 and HCN (jPonfria et al.l I2OO8I I 



to which their chernistry is closely coupled 
fsee lMillar k HerbstlEoM ) . C7H and CsH are 
also more abundant result of the addi- 
tion of CeH" to the chemical network which 
undergoes associative detachment with H to 
form C6H2. C6H2 reacts with C to produce 
C7H or with C2H to produce C8H2, which is 
photodissociated to CgH. 

Among the 33 chemical species listed in 
Table O 12 show calculated column densi- 
ties that differ from observations by more 
than a factor of 10. In the present model, 
the moderate-to-large-sized hydrocarbons and 
cyanopolyynes C4H2, C5N, CeH, CgH, and 
HCgN might be considered to have column 
densities that are discrepantly greater than 
the observations. By reducing the gas-to- 
dust ratio G to 50 (from the present value 
of 200), these discrepancies are eliminated 
due to the retardation of the C2H2 and HCN 



photochemistry as a consequence of increased 
dust extinction. The only species to be detri- 
mentally affected by such a reduction in G 
are C3N and C3N~, and as mentioned previ- 
ously, the chemistry involving these species is 
not completely understood at present. Tak- 
ing G = 50 may be too low however, be- 
cause values_^roin__tlieliter^ (see Table 
5 of iMen'shchikov et al.l boOlh are typically 
in the range 200-1000. It is plausible that 
other factors might contribute towards a re- 
duction in the photochemical yields, such as 
a re duced incident radiat ion field strength 
(see iMillar &: Herbsti 11994 ) and / or lower ini- 
tial C2H2 and HCN abundances, to bring the 
model results into better agreement with ob- 
servations. Changes in the rates of the reac- 
tions involved in the synthesis of moderate-to- 
large-sized hydrocarbons and cyanopolyynes 
could also help to resolve some of the discrep- 
ancies between the current model results and 
observations. For example, the reactions of 
the kind C2H + HC2„+iN HC2n+3N + H 
and C2H + C„H2 — > + H are crucial 

for the synthesis of successively larger species 
in the model, and the rates of these reactions 
are only known, at best, to within a factor 
of 2. Similarly, the photodissociation rates of 
many of the larger species in the model are 
calculated only from estimated cross-sections. 

It has been brought to our attention by the 
referee that in a recent study of narrow sub- 
millimetre emissio n lines from IRC +102 16 by 
Patel et al. ( 20081 ). a lower limit to the abun- 



dance of CS in the inner envelope of 9.3x 10~^ 
was derived. This is significantly greater than 
our adopted abundance of 4.0x 10~^. Thus we 
have investigated the impact on our model of 
raising the initial CS abundance to 1.0x10"^. 
We find that the only species to be apprecia- 
bly affected are those containing sulphur and 
those whose chemistries are closely related to 
the sulphur-containing molecules. The col- 
umn densities of HCS, OCS, H2CS and C„S 



16 



(for n = 1 — 5) scale approximately linearly 
with the initial CS abundance. The C3 and 
C3H column densities are raised by about 
50%, mainly result of the photodissoci- 
ation of C4S to produce C3 + CS. The OH 
column density increases by a similar fraction 
due to the hydrogen-transfer reaction of HCS 
with O. No other species in Table [3] have col- 
umn densitities that are significantly affected 
by the increase in the initial CS abundance. 

The inclusion of density-enhanced shells in 
the model has a profound impact on the cal- 
culated radial abundance distributions, both 
through the extinction effect of the shells 
which modifies the photochemical reaction 
rates, and through the effects of their in- 
creased densities which accelerate the binary 
chemical reactions. However, the impact of 
the shells on the total column densities is gen- 
erally rather small; even though the chemical 
abundances may vary by over an order of mag- 
nitude inside the shells, the fact that the shells 
are narrow compared to the radial extent 
of the CSE means that they contribute to- 
wards only a modest fraction of the total col- 
umn densities. Species whose chemistries are 
closely linked to the radiation field strength 
are most affected by the shells, for example 
C3H4 and CH3CHCN (whose dominant de- 
struction channels are by photodissociation), 
are shielded from dissociating radiation out 
to a larger radius by the increased extinction. 
The enhanced shielding in the envelope also 
reduces the ionisation rate, resulting in up to 
an order of magnitude reduction in the elec- 
tron and C"*" column densities at certain radii. 

A primary motivation of this study is to 
examine whether the struct ure observed by 
Dinh-V-Trung fc LimI (|2008l l in HC3N and 
HC5N maps of IRC-l-10216 is consistent with 
their suggestion that the gas and dust are cou- 
pled, sharing a similar density distribution. 
The inclusion of density-enhanced shells in 
our model results in peaks in the radial num- 



ber density distributions of molecules within 
the shells, which shows that density enhance- 
ments in the molecular gas are able to produce 
small-scale structure in the molecular distri- 
butions similar to those observed. 

When molecular excitation is taken into 
account, the emission intensity profiles for 
C2H, C4H, CeH and HC3N all peak inside the 
dense shell at 15". This finding is consistent 
wit h the 3 mm niaps of the se species made 
bv ICuelin et al.l (|l999l 'l and iLucas &: GuelinI 
(|l999l ) and shows that the presence of dense 
shells can cause the modeled molecular abun- 
dances to peak around the same radius and in 
a narrow region corresponding to the radius 
and thickness of the density-enhancement. 
The simulated emission maps (Figure [6]) show 
a similar picture for C2H and CgH. However, 
the inner broad intensity maximum of the 
C4H radial profile results in a ~ 5" wide, 
intense emission ring inside of the density- 
enhanced 15 " ring, which is not seen in the 



Guelin et al.l (|1999l ) map. Possible explana- 



tions for this discrepancy might be that the 
excitation of C4H favours its observation in 
the narrow shell (contrary to the crude ex- 
citation analysis we performed), or that the 
interferometer used in the observations failed 
to detect the broader structure of the inner 
ring. 

In the model, the gas and dust share 
identically-shaped density profiles which re- 
quires that the gas and dust must be dy- 
namically coupled . This is contr ary to 
the theory (see 



Truong-Bach et al. 1991 



Mauron &: Huggins that radiation pres- 

sure accelerates the circumstellar dust to a 
radial drift velocity ~ 2 km s~^ faster than 
the gas. Possible coupling mechanisms may 
include thermal and turbulent motion of the 
gas and dust. Future models may need to 
include a detailed analysis of the coupling be- 
tween the gas and dust shells. 

The shells used in the model represent a 
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gross simplification of the dust sliell structure 
observ ed in IRC+10216 by Mauron fc Hugging 
!3owever, the results presented here — 
that the maximum abundances of species con- 
taining carbon chains match the peaks in the 
circumstellar dust distribution — should be 
generally applicable to more complex repre- 
sentations of the density structure of the CSE. 
The density-dependent nature of molecular 
excitation means that detailed excitation cal- 
culations are required to determine the rela- 
tionship between observed emission intensity 
and molecular abundance distributions. In 
our model the temperature distribution of the 
gas is assumed to be continuous, with the as- 
sumption that the density-enhancements are 
in thermal equilibrium with the surrounding 
medium. More detailed observations of the 
gas and dust will be required to determine if 
this assumption is realistic. In addition, ac- 
curate collisional and vibrational excitation 
rates will be required in order to more accu- 
rately calculate the rotational excitation of 
the molecules considered. 

5. Conclusion 



Density-enhanced shells of gas and dust 
have a significant impact on the calculated ra- 
dial distributions of molecules in the expand- 
ing envelop e of IRC+10216. Based on th e sug- 
gestion by IPinh-V-Trung k LimI (|2008l ) that 
the gas and dust are coupled, we included in 
a new model for the CSE a set of density- 
enhancements with para nieters based on the 
dust s hell observations bv lMauron Sz Huggins 
Photochemistry is delayed to greater 
radii due to increased shielding by the dust 
shells. Density-enhancements in the gas re- 
sult in molecular abundances that peak inside 
the narrow shells, showing that the clumpy 
structure observed in IIC3N and IIC5N may 
be the result of density enhancements in the 
molecular gas. The calculated emission in- 
tensity profiles for C2H, CgH and IIC3N all 



peak within a narrow band about 15" from 
the central star, which is consistent with de- 
tailed emission maps of these species. The 
emission profile for C4H has a broad maxi- 
mum around 8" that is not present in observed 
maps, which may be indicative of a need for 
further improvements in the model. 

The ethynyl anion C2H~ and the nitrile an- 
ions C2n-iN~ (for n = 1,2,3,4), have been 
calculated to reach observable abundances in 
the CSE. For the smaller nitrile anions for 
which radiative electron attachment is slow 
(i.e. CN~ and CaN"), this result is uncertain 
due to the unknown product branching ratios 
in the associative electron detachment chan- 
nels of the reactions N + C~ and N + CmH~. 
Further laboratory measurements (including 
up to at least m = 10), will be required in 
order to confirm the importance of these re- 
actions in anion astrochemistry. C2H~ is pre- 
dicted to be produced in abundance in the in- 
ner CSE (at much smaller radii than the other 
anions in the model), product of the reac- 
tion of H" with C2H2. Other possible proton 
transfer reactions of this kind {i.e. XH + Y~ 
— > X~ + YH, for the anions listed in Ta- 
ble [2]), also need to be considered for possible 
inclusion in future anion chemical networks. 

Our chemical models produce MgNC with 
a peak abundance in the outer envelope. De- 
pendent on the initial Mg abundance used, the 
observed MgNC column density matches ob- 
servation, which shows that gas-phase chem- 
istry is a viable route to the formation of this 
species. 

We gratefully acknowledge Veronica Bier- 
baum and Eric Herbst for their contribu- 
tions to this work regarding the calculation 
of branching ratios for reactions between an- 
ions and nitrogen atoms. MAC thanks QUB 
for financial support. Astrophysics at QUB is 
supported by a grant from the STFC. 
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